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Aims. We investigate the performance of an analytic model of the 3D matter distribution, which combines perturbation 
theory with halo models, for weak-lensing configuration-space statistics. 

Methods. We compared our predictions for the weak-lensing convergence two-point and three-point correlation functions 
with numerical simulations and fitting formulas proposed in previous works. We also considered the second- and third- 
order moments of the smoothed convergence and of the aperture-mass. 

Results. As in our previous study of Fourier-space weak-lensing statistics, we find that our model agrees better with 
simulations than previously published fitting formulas. Moreover, we recover the dependence on cosmology of these 
weak-lensing statistics and we can describe multi-scale moments. This approach allows us to obtain the quantitative 
relationship between these integrated weak-lensing statistics and the various contributions to the underlying 3D density 
fluctuations, decomposed over perturbative, two-halo, or one-halo terms. 
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1. Introduction 

The standard paradigm known as the ACDM cosmology 
includes dark compo nents: dark matter and dark energy 
(|Komatsu et al.|[2011[ ). Weak lensing of background galax- 
ies by foreground large-scale structures, the so-called "cos- 
mic shear" , has been recognized as a potentially pow- 
erful tool for probing the distribution o f dark matter 
as w ell as the nature of dark energy (jAlbrecht et al.l 
I2006D . Reports of significant detections of cosmic 
shear have been made by vario us groups (iBacon et al.l 



Hamana et al. 



Van Waerbeke et al 
200l 



200a iFu et al. 



200C; 



Jarvis et al 



[2008FlSchrabback et al.ll201ol) 



Wittinan et al.l l2000t 
20061 ISemboloni et al.l 



By analyzing the cosmic shear data, one can di- 
rectly measure the power spectrum of the matter den- 
sity fluctuations on cosmological scales, which con- 
tain a wealth of cosmological information such as neu- 
trino masses and dark energy equation-of-state parame- 
ters (Jarvis et al."2006; 'Semboloni et al.l 120061: llchiki et al.l 
12009; .Schrabback et a l. 2010). Consequently, it IS a mam 
goal of cosmology to infer and constrain these quanti- 
ties from observations. To do this, a number of ambi- 
tious surveys are planned, s uch as the Hyper Su prime- 
Cam Weak Lensing Survey (jMivazaki et al.l l2006l n. the 
Dark Energy Survey (DESj3, the Large Synoptic Survey 



Telescope (LSSTfl Euchd (iRefregier et al.ll2nifllfl and the 
Wide-Field Infrared Survey Telescope (WFIRStB- 

Most weak-lensing information is contained in small an- 
gular scales and therefo re weak-lensing stat i stics are nonlin- 
ear and non-Gaussian (iMunshi et al.ll2008l : iTakada fc JainI 
[20091 : [Satoet al.ll2009L l2011bD . If we aim to exploit the full 
information, we have to treat the nonlinear effects to ac- 
curately model the weak lensing statistics. Furthermore, 
one has to use an appropriate likelihood function with 
given marginal distributions, oth erwise the obtained resu lts 
would be systematically biased (jSato et al.l |2010l l2011aD. 

In a first companion paper ( Valagc as et al.l l2011bL 
hereafter paper I), we studied the Fourier-space weak- 
lensing statistics such as the weak-lensing power spectrum 
and bispectrum, and found th at our model proposed by 
IValageas &: Nishimichil ()2011al lbl). which combines pertur- 
bation theory with halo models, agrees better with ray- 
tracing simulations than previously published fitting for- 
mulas and phenomenological models. In this second paper, 
we study the real-space weak-lensing statistics, which are 
more often used for the statistical analysis of actual mea- 
surements than Fourier-space statistics, because observa- 
tions are made in configuration space. 

Previous works have already shown that on small 
scales the halo model provides a good description of the 
two-, three- and four-point correlations or smoothed mo- 
ments of the cosmic shear (using some approximations ) 
(jTakada fc Jainll200ll2003t iBenabed fc Scoccima^l2006t ). 
whereas a stochastic halo model can recover the prob- 



^ http:/ /www. naoj.org/Projects/HSC/index.html 
^ http:/ /www. darkenergysurvey.org/ 



^ http://www.lsst.org/ 

^ http://www.euclid-ec.org/ 

^ http://wfirst.gsfc.nasa.gov/ 
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ability distribution function of th e unsmoothed conver- 
gence ()Kainulainen fc Marrall2011allb[ ). Here we include all 
"one-halo", "two-halo" and "three-halo" terms, as well as 
one-loop perturbative results, and we compare these with 
larger-scale simulations. This yields a greater accuracy and 
allows us to compare these different contributions, from 
very large to small scales. This should be useful for prac- 
tical purposes because these different terms have different 
theoretical accuracies and probe different regimes of gravi- 
tational clustering, hence it is important to know their rela- 
tive impact on weak-lensing probes as a function of angular 
scale. 

This paper is organized as follows. In Sect. [2] we briefly 
recall how configuration-space weak-lensing statistics are 
computed from polyspectra of the 3D matter density field. 
We describe our numerical simulations and the data anal- 
ysis in Sect. [31 Then, we present detailed comparisons be- 
tween the simulation results, previous models, and our the- 
oretical predictions for two-point functions in Sect. |4] and 
three-point functions in Sect. [SJ We study the relative im- 
portance of the different contributions arising from "one- 
halo", "two- halo", or "three-halo" terms in Sect. [6l Then, 
we check the robustness of our model when we vary the cos- 
mological parameters in Sect. [7] and we briefly study multi- 
scale moments in Sect. [H Finally, we conclude in Sect. [9] 

2. From 3D statistics to weak-lensing statistics 

2.1. Lensing power spectrum and bispectrum 

Using Born's approximation, the weak-lensing convergence 
k{6) can be written as th e integral of the density con- 
trast along the line of sight ([Bartelmann fc Schneideijr2001t 
iMunshi et alll2008[ ). 



<0) 



(1) 



where x ^'^d 2? are the radial and angular comoving dis- 
tances, 

2c^V{xs) 

and Zs is the redshift of the source (in this article we only 
consider the case where all sources are located at a single 
redshift to simplify the comparisons with numerical simu- 
lations and the dependence on the source redshift). Then, 
using a flat-sky approxim ation, which is valid fo r small an- 
gles below a few degrees (jValaeeas et al.l[2bllal ). we define 
its 2D Fourier transform through 



(3) 



As in paper I, we define the 2D convergence power spectrum 
and bispectrum as 

{K{ei)R{£2)) = SD(£l+£2)PAil), (4) 

and 



{k{£i)h{£2)k{£3)) - SD{£l+£2+£3)BAeij2,h)- 



(5) 



From Eg .(IT]) one obtains at once f rom Limber's approxima- 
tion (Limber 1953; Kaiseij [l992l : iBartelmann fc Schneideij 
2001; Munshi et al... 200S ) 



P.{e) = 27r 



5.(^1,^2,^3) = (27r)2y^ dx:^B{h/V,e2/V,e3/V;z), (7) 

where P{k;z) and i?(fci, fc2, fca; 2) are the 3D power spec- 
trum and bispectrum of the matter density contrast at 
redshift z. As described in paper I, this provides the 
weak-lensing convergence power sp ectrum and bispectrum 
from t he model we developed in IValageas fc Nishimichil 
(l2011al lb[) for the 3D power spectrum and bispectrum 
through a simple integration over the radial coordinate up 
to the source plane. 



2.2. Configuration-space statistics 

We focus here on configuration-space weak-lensing statis- 
tics, which may be more convenient than Fourier-space 
quantities for practical purposes because of complex sur- 
vey geometries. Indeed, observations of the shear field are 
made in configuration space, by measuring large-scale cor- 
relations of galaxy ellipticities, and going to Fourier space 
(with further operations that are not local in real space) 
can be difficult because galaxy surveys do not cover the 
whole sky and shear maps can show irregular boundaries 
and internal holes due to observational constraints. 

In particular, we consider the two-point and three-point 
correlation functions of the weak-lensing convergence, de- 
fined as 

{K{e,)K{e2)) =u\d2-ei\), (8) 

(k(0i)k(02)aj(03)) - U\02 - e-sl |03 - 0l\, Wl - 02|), (9) 

where we used statistical homogeneity and isotropy (hence 
Ck only depends on the lengths of the three sides of the 
triangle defined by the summits {61,92,63} and the bis- 
pectrum d?]) only depends on the three lengths {^1,^2,^3})- 
These real-space correlations can be expressed in terms 
of the Fourier-space power spectrum and bispectrum as 



U0) = 2n d£eP^{£)M£9), 

JQ 



(10) 



where 9 = 1^2 ^ ^1 1 is the pair angular distance as in Eq.([i 
and 

00 /'^/2 /'27r 

d£f / d9 sm{2e) / d^ 

^1*^2 "'0 Jo Jq 



C,k{vi,V2,V3.) 



xB^{£i,£2,£3) Jo(Vl + sin(20) cos((^ + ag)), (11) 

where {j^i, !^2, 1^3} are the lengths of the three sides of the 
triangle {61,62,63} as in Eq.Q, as is the inner angle at 
summit 63, we chose without loss of generality 



'^i>i^2> V3, cosaa 



(12) 



dx^P{£/V;z), 



(6) 



2viV2 

and we noted in Eq. pTj) the multipoles 

= i£^, £2-^, £l^£\+£l + 2£,£2Cos^. (13) 

V2 Vl 

In addition to the three-point correlation C,^ of the 
convergence, it can be useful to consider the three-point 
correlation of the cosmic shear 7. Because the latter is 
a spin-2 quantity, one is led to consider several three- 
point shear correlations (or "natural components"), de- 
pending on the choice of the reference direction (or of 
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the projection procedure of the cosmic shear vectors) , 
see [Schneider fc Lomba^ (|2003D ^Schneider et all (|2005li . 
However, they can all be written as integrals over the con- 
vergence bispectrum, such as Ea. ljlip . or i ntegrals over th e 
convergence three-point correlation ^ (jShi et al.l 1201 If) . 
Therefore, although we only consider the convergence three- 
point correlation ([9|) in this paper, we can expect a similar 
level of agreement between our analytical model and simu- 
lations for these other three-point correlations (in addition 
we also consider the third-order moment of the aperture- 
mass, which can be related to both the convergence and the 
cosmic shear). 

It is also common practice to study smoothed averages 
Xs of the convergence or shear field, defined by their filter- 



ing window Wg ° (9) through 

/c\f) 
— n{e)w,f{e). 



(14) 



For instance, the smoothed convergence is defined by a 
top-hat filtering, 



Wg^iO) = 1 if \0\ < 9s and zero otherwise, 



(15) 



while the "aperture-mass" Map is defined by a com pensated 
filter (ISchneideilllQQilVan Waerbeke et alJl200l[ ). such as 



W^'^-^iO) =3(1 



1-3- ) if \e\<es, (16) 



and lYg'^'^^ie) = if |0| > 61,. This also reads in Fourier 
space as 



with 



de 



In particular, we have 
Ji 



s) =2- 



(17) 



(18) 



(19) 



We mainly focus here on one-point moments of and 
Map, that is (XP), and we do not consider multi-point 
statistics such as {Xs{9i;9si)..Xs{9p; 6sp)) associated with 
p windows centered on p different directions and with p dif- 
ferent angular radii. However, we will check the validity of 
our model for multi-scale statistics, that is, for windows of 
different sizes centered on the same direction, in Sect. |8l 

In the simpler case of one-point moments (X^) the vari- 
ance reads as 



Jo 



(20) 



while the third-order moment reads as 



poo nil /"arccosf^i/ {2^2)] 

247r / d£i£i / d^2^2 / dip 

Jo Jii/2 Jo 

xB«(^i, ^2, 4) w^; iWWg''; {i20s)w^; (W. (21) 

where we used the symmetries of the bispectrum and we 
noted 



el = el + ij ~ 2ii£2 cos If. 



(22) 



In practice, to avoid the numerous oscillations and changes 
of sign brought by the Fourier-space filters W^' given in 
Ea.(|19p. we found it convenient to express the third-order 
moment (1211) in terms of the real-space three-point corre- 
lation (|9]), although this yields a five-dimensional integral 
instead of the three-dimensional integral (I^TI) . 



{X's 



with 





2tt 



X / de^esWg'l^ies) / d(^2 / d^3CK(^i,^2,j^3),(23) 

Jo ° Jo Jo 



93 COS (^3 - 6*2 COS 1^2)^ + (03 sin 1^3 - 62 sin 932)^, (24) 



^3 cos ip3 
92 COS ip2 



i)^+0i sin2 (^3, 



9l sin^ ip2. 



(25) 
(26) 



3. Numerical simulations 



We performed the ray-tracing simulations through high- 
resolution A^-body simulations of cosmological structure 



formation ("Jain et al. 


2000l:lHamana et al.ll200lUSato et al.l 


2009: Takahashi et al. 


|201l|) to obtain accurate predictions 



of the configuration statistics for weak lensing. To run the 
iV-body simulat ions, we used a modified version of the 
Gadget-2 code (jSpringell 120051 ) and employed 256^ parti- 
cles for each simulation. The ray-tracing simulations were 
constructed from 2 x 200 realizations of A^-body simula- 
tions with cubic 240 and 480/i~^Mpc on a side, respec- 
tively, to cover a ligh t cone of angular size 5° x 5° (see 
Fig. 1 in ISato et al.l l2009). For our fiducial cosmology, we 
adopted the standard ACDM cosmology with matter frac- 
tion fim = 0.238, baryon fraction fib = 0.0416, dark energy 
fraction Ode = 0.762 with the equation of state parame- 
ters Wo = —1 and Wa = 0, spectral index Ug = 0.958, 
normalization Ag — 2.35 x 10~^, and Hubble parameter 
h = 0.732, which are consistent with the WMAP three- 
year results (IS oergel et al.l [20071 ). This fiducial cosmology 
gives the normalization erg = 0.759 for the variance of the 
linear density fluctuations in a sphere of radius 8/i~^Mpc. 
We considered source redshifts at either Zs ~ 0.6, 1.0, or 
1.5. Using ray-tracing simulations we generated 1000 real- 
izations of convergence maps for each source redshift. 

In addition to the fiducial cosmology case, we performed 
ray-tracing simulations for several slightly different cos- 
mologies. We varied each of the following cosmological pa- 
rameters: As, Us, the cold dark matter density Q,ch^, ^dc, 
and Wo by ±10%, respectively. Therefore h, r2,„, and fib 
are dependent parameters, because we assumed that the 
Universe is fiat and the baryon density fib^-^ is fixed. For 
each of these ten different cosmologies, we obtain 40 real- 
izations of convergence fields for each of the three source 
redshifts. Details of the me thods used for th e ray-tracing 
simulations can be found in ISato et al.l ([20091) . All conver- 
gence maps used in this paper are the same as those used 
in paper I. 

In Sects. 01111 we use the maps for the fiducial cosmology, 
while in Sect. [71 we use those for the varied cosmologies 
to investigate the robustness of our model. In Sect. |7l we 
show the results for six cases, varying n,, ftch^, and wq by 
±10%. The exact values for these cosmological parameters 
are listed in Table A.l in paper I. 
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Fig. 1. Upper row: convergence two-point correlation function for sources at redshifts Zs = 0.6, 1, and 1.5, as a function 
of the angular pair separation 6. The points are the results from numerical simulations with 3 — cr error bars. The low 
bl ack dashed line "L" is the linear correlation, the middle blue dash-dotted line "S" is the result from the "halo-fit" 
of ISmith et al.l (|2003l ). and the upper red solid line "comb." is the result from our model, which combines one-loop 
perturbation theory with a halo model. The vertical arrow shows the scale down to which the simulation result is valid 
within 5%. Middle row: Variance of the smoothed convergence for the same cases, as a function of the smoothing angle 
9s. Lower row: Variance of the aperture mass for the same cases, as a function of the smoothing angle 9^- 



4. Lensing two-point functions 

We now compare our results for weak-lensing two-point 
functions with numerical simulations. As in paper I, we 
also considered the predictions obtained from the popular 
"halo-fit" fitti ng function for the 3D density power spec- 
trum given in ISmith et al .' ('2003'), to estimate the advan- 
tages of more systematic approaches like ours. 

We show our results for the convergence two-point cor- 
relation £,K.{9), the variance of the smoothed convergence 
{k^}, and the variance of the aperture mass (M^^) in Fig.[T] 
The numerical error bars increase on large scales because of 



the finite size of the simulation box. On small scales the nu- 
merical error is dominated by systematic effects, because of 
the finite resolution, which leads to an underestimation of 
the small-scale power. This underestimation was clearly ap- 
parent for the power spectrum Pk(^) shown in paper I and 
can also be seen (especially at low redshift) for the variance 
of the aperture mass, which involves a compensated filter 
and probes a narrow range of scales (|Schneidedll"996l) . For 
each source redshift we estimated the angular scale down 
to which the simulations have an accuracy of better than 
5% by comparing with higher-resolution simulations (with 
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512^ particles instead of 256"^). This scale is shown by the 
vertical arrow in Fig. [T] and we can check that our model 
indeed agrees with the numerical simulations down to this 
angular scale. 

The two-point correlation and the smoothed conver- 
gence are not as sensitive to this low-resolution effect be- 
cause they involve uncompensated filters, W{0) ^ 0, which 
implies that at a given smoothing angular scale they 
receive greater contributions from larger scales (which are 
unaffected by the numerical resolution) than the aperture 
mass. 

For the same reason, ^^(0) and (k^) remain adequately 
described by linear theory down to 10 arcmin, whereas 
(Mfp) already shows significant deviations at ^ 100 ar- 
cmin. This also explains why the predictions from our 
model and the "halo-fit" are closer for the first two statis- 
tics than for Afap. In agreement with our results for the 
convergence power spectrum (paper I) and previous works 
(iWhite fc VaS[200llHilbert et al.l [2009t iSato et al.l[2009t 
iSemboloni et al.r 20111) ■ the "halo- fit" formula underesti- 



mates the power on small scales. The discrepancy is again 
larger for the aperture mass. As in paper I, our model 
provides a more accurate match to the numerical simu- 
lations down to ~ 1 arcmin. On smaller scales the re- 
sults from simulations show a fast drop (especially for 
Map) that is not physical but due to the finite resolu- 
tion. This prevents an accurate comparison with our pre- 
dictions. However, since our model is built from a physical 
halo model and has been tested for the 3D density field 
down to highly nonlinear scales with hi gher-resolution sim- 
ulations H Valageas fc Nishimichil [20 llaf) . it should be more 
reliable than the numerical results shown in Fig. [T] below 
~ 1 arcmin. This shows one advantage of analytic (or semi- 
analytic) approaches as compared with numerical simula- 
tions: they can provide realistic predictions on a wider range 
of scales. 



5. Lensing three-point functions 

We now compare our analytical results with numerical 
simulations for three-point functions. As in paper I, we 
also considered the predictions obtained from the following 
three simple models, which have been used in some previous 
works. ( We did not consider the scale transformation intro- 
duced in iPan et all (|2007i) because we have already shown 
in paper I that it does not provide a sufficiently accurate 
model for the convergence bispectrum. We have checked 
that we obtain similar results for (k^) and (M^p).) 

The first model, "treet", is the lowest-order ("tree- 
order") prediction from standar d perturbation theory , 
which reads for the 3D bispectrum (jBernardeau et aO2002[ ) 

5trcoJfci,fc2,fc3)=2F2(fci,fc2,m2)PL(fcl)^'L(fc2) + 2cyC.(27) 



where [in = (ki • k2)/(fcifc2) and 



-^2(^1,^2,^12) 



1 / fci ki 



2 V ^2 fci 



M12 



■M?2- 



(28) 



The second model, "trecNL", is given by Eq. ipT)) where we 
replace the linear 3D power P L(k) by the nonlinear power 
Ps{k) from lSmith eTall (l2003ll . The third model, "F2,nl", 
makes the additional modification to replace the kernel F2 
by an effective kernel -F2.NL that interpolates from the large- 
scale perturbative result ((28)) to a small-scale ansatz where 



the angular d ependence on /j,i2 vanishes , using the fitting 
formula from IScoccimarro fc CouchmanI (|2001 ). We con- 
sidered two variants, using ei ther the "halo-fi t" nonlinear 
power spectrum Ps{k) from Smith et al.l ((2003) or the non- 
line ar power spectrum Ptang(fc) of our model, see paper I 
and lValageas fc Nishimichil ( 2011bl ). 

We show our results for the convergence three-point cor- 
relation function (for equilateral triangles) and for the 
third-order moments (k^) and (M^p) in Fig. [5| In agree- 
ment with the second-order statistics shown in Fig. [1] the 
lowest-order perturbation theory prediction, "trecL" , re- 
mains valid down to smaller angular scales for Ck and Ks 
than for Map because the former involve uncompensated 
filters instead of a compensated filter. 

The simulations slightly underestimate the power on 
moderate and large angular scales because of the finite size 
of the simulation box (240/i~^ and 480h~^ Mpc, which cor- 
responds to 343 and 686 arcmin at z — 1), which cuts con- 
tributions from longer wavelengths. This discrepancy is not 
due to the analytic model because on the largest angular 
scales we can check that all theoretical predictions converge 
on the linear theory (as they must for CDM power spec- 
tra) while predicting somewhat more power than measured 
in the simulations (see the first row in Fig. [2|). Therefore, 
this mismatch is not caused by higher-order perturbative 
corrections (e.g., two-loop terms), which are even smaller 
than the one-loop contributions that we included in our 
model. In agreement with this explanation, the discrep- 
ancy is smaller for Afap than for and k because, for a 
given smoothing radius 9s, Afap is less sensitive to larger 

scales thanks to its compensated filter Wg^'^^ . This shows 
that analytical models, such as the one we propose here, 
are competitive with numerical simulations if one needs to 
describe a broad range of scales. 

In agreement with the results obtained in paper I for the 
convergence bispectrum, using the nonlinear power within 
the "tree-level" expression ^17} significantly increases the 
third-order moments on smaller scales and improves the 
general shape but is not sufficient to bridge the gap with 
the simulations . Modifying the kernel F2 by usi ng the fit- 
ting formula of IScoccimarro fc Couchma^ (l2001h improves 
the predictions even more, especially when we use the 3D 
nonlinear power spectrum given by our model, which was 
shown ear lier to be reasonably accurate (s ee Fig. [1] and pa- 
per I and IValageas fc Nishimichil (|2011bD ). I ndeed, as for 
the va riance, the "halo-fit" power spectrum of ISmith et al.l 
(|2003l ) yields too littl e power on small scale s, in agreement 
with previous works ([Semboloni et al.ll2011l) . However, this 
approach still underestimates the weak-lensing signal. 

The best agreement with the numerical simulations is 
provided by our model. As was seen for the convergence 
bispectrum in paper I, it is interesting to note again the 
good match on the transition scales, ~ 5 arcmin, which 
a priori are the most difficult to reproduce since they are at 
the limit of validity of both perturbative approaches (which 
break down at shell crossing) and halo models (which as- 
sume virialized halos). On small angular scales we again 
predict more power than is measured in the simulations, 
but like for the second-order moments this is at least partly 
caused by the lack of small-scale power in the simulations 
because of the finite resolution. Thus, we again plot in 
Fig- m the vertical arrows that were plotted in Fig. [TJ Since 
the bispectrum typically scales as the square of the power 
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Fig. 2. Upper row: convergence three-point correlation function for equilateral triangles, as a function of the triangle 
side u, for sources at redshifts Zs = 0.6, 1, and 1.5. The points are the results from numerical simulations with 3 — cr 
error bars. The low black dashed line "treeL" , the black dotted line "treeNL" and the two blue dashed lines "F2_nl" are 
obtained from the ansatz (|27p . using either the linear 3D power spectrum, or the "halo-fit" power, or an effective kernel 
-Pj.NL with the "halo-fit" power (lower curve) or the power from our model (upper curve). The red solid line "comb." 
is our combined model, described in paper I. The vertical arrows are at the same angular scale as in Fig. [TJ Middle 
row: Third-order moment of the smoothed convergence, as a function of the smoothing angle 9s . Lower row: Third-order 
moment of the aperture mass, as a function of the smoothing angle 9s. 



spectrum, this should roughly correspond to an accuracy 
threshold of about 10% for the simulations. We can check 
that our model agrees with the numerical results down to 
this angular scale. As in Fig.[Tl Map is much more sensitive 
than and Kg to this finite-resolution effect. 



6. Relative importance of the different 
contributions 

As in paper I, we have seen in the previous section that our 
model, which is based on a combination of perturbation 
theories and halo models, provides a good match to numer- 
ical simulations. Therefore, it can be used to predict weak- 
lensing statistics for a variety of cosmologies, which is an 
important goal for observational and practical purposes. A 
second use of our approach is to compare the different con- 
tributions that eventually add up to the signal that can be 
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Fig. 3. Convergence two-point correlation function (left panel), smoothed convergence (middle panel), and aperture- 
mass (right panel) for sources at redshift Zs = 1. The points are the results from numerical simulations with 3 — ct error 
bars. The low black dashed line "L" is the linear prediction, the middle black dash-dotted line "1-loop" is the two-halo 
contribution, for which we used a perturbative resummation that is complete up to one-loop order, the upper blue dashed 
line "IH" is the one-halo contribution, and the red solid line is our full model, as in Fig. [TJ 




Fig. 4. Convergence three-point correlation for equilateral configurations (left panel) and third-order moments (k^) 
(middle panel) and (Mfp) (right panel), for sources at redshift Zg = 1- The points are the results from numerical 
simulations with 3 — cr error bars. The low black dashed line "L" is the lowest-order perturbative prediction (|27p . the 
"l-loop" dash-dotted black line is the prediction of one-loop standard perturbation theory, which is identified with our 
three-halo term, the blue dotted line "2H" is the two-halo contribution and the upper blue dashed line "IH" is the 
one-halo contribution. The red solid line is our full model. 



measured in weak-lensing surveys. Thus, we can distinguish 
the various perturbative terms as well as the nonperturba- 
tive contributions associated for instance with one-halo or 
two-halo terms. This is a second advantage of these mod- 
els, as compared with fitting formulas or direct ray- tracing 
simulations. This enables a deeper understanding of which 
properties of the matter distribution, and of the overall cos- 
mological setting, can be probed by a given gravitational 
lensing measure. It can also help to estimate the accuracy 
that can be aimed at in weak-lensing statistics as a func- 
tion of scales, because different contributions suffer from 
different theoretical uncertainties. 

6.1. Two-point statistics 

We plot our results for (k^), and (Mf ) in Fig.[3]at red- 
shift Zs = 1. In addition to the full model prediction that 
was already shown in Fig.[l] we show the underlying 2-halo 
and 1-halo contributions. Because of their uncompensated 



filters, which make statistics at a given angular scale re- 
ceive contributions from 3D fluctuations on a wide range of 
scales, taking the one-loop perturbative term into account 
only yields a small increase of and (k^) over a wide range 
of angular scales, as compared with the linear prediction. 
In contrast, for (Mfp) this one-loop contribution peaks on 
a narrow range of angular scales around 40 arcmin and has 
a significant impact that improves the match to the nu- 
merical results. On small scales the two-halo contribution 
decreases close to the linear prediction thanks to the partial 
resummation of hi gher perturbative orders. As explained in 
IValageas fc Nishimich i (2 011a) and paper I, this is a useful 
improvement over the standard one-loop perturbation the- 
ory because it ensures that the two-halo term does not give 
significant contributions on very small scales, in agreement 
with physical expectations. 

Then, these two-point weak-lensing quantities become 
dominated on small scales by the one-halo term but like for 
the 3D and 2D power spectra, there remains a significant in- 
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termediate range. Again, our model provides a satisfactory 
interpolation on these scales, but it would be interesting to 
build a more systematic procedure, for instance by includ- 
ing higher orders of perturbation theory or by building a 
more refined matching between the two-halo and one-halo 
regimes. In any case. Fig. [3] clearly shows how (Ks), and 
(M^p) depend on large-scale perturbative density fluctua- 
tions or on small-scale halo properties, as the angular scale 
varies. 

We checked that we obtain similar results at redshifts 
Zs = 0.6 and 1.5. 

6.2. Three-point statistics 

We now study the various contributions to the lensing 
three-point functions, associated with the thee-halo, two- 
halo, and one-halo terms. 

We plot our results for (for equilateral configura- 
tions), (k^), and (Mfp) in Fig. HJ at redshift Zs = 1. As 
in paper I, the three-halo term is identified with the per- 
turbative prediction and we used the standard perturba- 
tion theory at one-loop order. Like for the 3D bispectrum 
(| Valageas fc Nishimichil[2011b[ ). and contrary to the power 
spectrum, this gives a contribution that becomes negligi- 
ble on small scales, so that it is not necessary to use a 
resummation scheme or to add a non-perturbative cutoff 
to ensure a good small-scale behavior. On the other hand, 
contrary to the two-point statistics shown in Fig. [31 going to 
one-loop order now provides a great improvement over the 
tree-order result, even for and (k^) . This feature was al- 
ready noticed for the 30 bisp ectrum (jSefusatti et al.l[20Tot 
I Valageas fc Nishimichill20TTbl ) and the convergence bispec- 
trum (paper I). Thus, combining this one- loop perturba- 
tive contribution with the two-halo and one-halo terms is 
sufficient to obtain a good match to the simulations, from 
the quasilinear to the highly nonlinear scales. This suggests 
that higher orders of perturbation theory do not signifi- 
cantly contribute to the bispectrum and that we already 
have a reasonably successful model. The two-halo term is 
also subdominant on all scales (by a factor ^ 10 at least). 
This is a nice property since because it is a mixed term, 
which involves both large-scale halo correlations and in- 
ternal halo structures, it may be more difficult to predict 
than the three-halo term (which is derived from system- 
atic perturbation theories) and the one-halo term (which 
only depends on internal halo profiles and mass function). 
These variou s features were also observed for the 3D bispec- 
trum (jValag eas & Nishimichi 2011b) and the convergence 
bispectrum (paper I). 

Again, we checked that we obtain similar results at red- 
shifts Zs — 0.6 and 1.5. 



7. Dependence on cosmology 

In this section we check the robustness of our model when 
we vary the cosmological parameters. As in paper I, we con- 
sidered six alternative cosmologies, where ris, flch^, and wq 
are modified by ±10% with respect to the fiducial cosmol- 
ogy used in the previous sections. The values of the associ- 
ated cosmological parameters are given in Table I in App. A 
of paper I. We compare the predictions of our model with 
numerical simulations for these six alternative cosmologies 
in Figs. [5] and [6] for the two-point and three-point statistics 



at Zs = 1 . To avoid overcrowding the figures we did not plot 
the error bars of the numerical simulations. Each pair n^, 
rich'^, and wq gives two curves that are roughly symmetric 
around the fiducial cosmology result, because we consider 
small deviations of ±10%. The deviations are largest for 
the ris case, which changes the shape of the initial power 
spectrum as well as the normalization erg. These six cases 
roughly cover the range that is allowed by current data, and 
the Us c ase is already somewh at beyond the observational 
bounds (jKomatsu et al.l [20111 ). Therefore, they provide a 
good check of the robustness of our model for realistic sce- 
narios. 

Like for the Fourier-space statistics studied in paper I, 
the dependence on cosmology of the two-point statistics is 
well reproduced by our model. For the three-point statis- 
tics it is not easy to make a very precise comparison be- 
cause the numerical results show a greater level of noise 
and are sensitive to finite resolution and finite size effects. 
However, where the simulations are reliable, we also obtain 
a good match with our predictions. We obtained similar 
results for Zg = 0.6 and Zg = 1.5, as well as for other cos- 
mologies where we vary Ag or fide by ±10%. This shows 
that our model and, more generally, models based on com- 
binations of perturbation theory and halo models provide 
a good modeling of the matter distribution and of weak 
gravitational lensing effects and capture their dependence 
on cosmology. Moreover, Figs.[5]and|n]clearly show that this 
analytical modeling is competitive with current ray-tracing 
simulations, because it provides reliable predictions over a 
greater range of scales. In particular. Figs. [5] and [5] show 
that the accuracy of our model is sufficient to constrain n^, 
^ch^, and wq to better than 10%. 

8. Multi-scale moments 

In the previous sections we considered single-scale mo- 
ments, {Xs{9sy), associated with one smoothing window 
with a single angular radius Og. One can also use 
multi-point statistics such as {Xs{0i;9si)--Xs{6p;Osp)) as- 
sociated with p windows centered on p different directions 
6i and with p different radii 9si ■ In this section, we briefly 
check the validity of our model for the case of centered 
multi-scale moments. To do this, all quantities X^i are cen- 
tered on the same direction on the sky but we allow the 
angular radii 9 si to be different. Then, Eas. (P0| and (1^31) 
generalize as 

/>oo 

{Xs{9si)Xs{9s2))^2tt / diep^ie) w^;{e9si)w^;{£9s2) 

jQ 

(29) 

and 

{Xs{9si)Xs{9s2)Xs{9si)) - 47r / -^W^l{9i) 

X / d(^2 / Aip3,(^:^{vi,V2,V3,). (30) 

We show in Fig. [7] the two-scale second-order moments 
{Ks{9s)K.s{a9s)) and {M^p{9s)M^^{a9s)) for a scale-ratio 
a = 2,5, and 10, at Zs — 1. A higher a yields a smaller 
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Fig. 5. Convergence two-point correlation function (left panel), smoothed convergence (middle panel), and aperture- mass 
(right panel) for sources at redshift — 1 for six cosmologies. The points are the results from numerical simulations and 
the lines are the predictions of our model. 




Fig. 6. Convergence three-point correlation for equilateral configurations (left panel) and third-order moments (k^) 
(middle panel) and (M^^) (right panel) for sources at redshift Zg = 1 for six cosmologies. The points are the results from 
numerical simulations and the lines are the predictions of our model. 



moment because it corresponds to a larger second angu- 
lar radius a9s. We obtain the same level of agreement as 
for the single-scale variances shown in Fig. [T] In particu- 
lar, we obtain a good match on small angular scales where 
the "halo-fit" formula somewhat underestimates the weak- 
lensing power. On large scales our analytical results are 
somewhat larger than the data obtained from the numerical 
simulations. This is due to the missing of large-scale modes 
in the simulations because of the finite size of the simula- 
tion boxes (240/i~^ and 480/i~^ Mpc) that correspond to 
343 and 686 arcmin at 2; = 1). This is more apparent than in 
the single-scale plots of Fig.[T]because we probe larger scales 
since the factor a is greater than unity. This again shows 
the advantage of analytical models such as ours, which are 
competitive with current ray-tracing numerical simulations 
to describe a broad range of scales. 

We show in Fig. [5] the three-scale third-order moments 
{Ks{9s)Ks{ae,)Ks{l3es)) and (Map(0.)Afap(a0.)Map(/30.)) 
for the scale-ratios {«,/?} — {2,5}, {3,9}, and {5,10}, at 

= 1. Higher values of a and /? give a smaller moment 
since they correspond to larger second and third angular 
radii a9s and /39s- We obtain the same level of agreement 
as for the single-scale third-order moments shown in Fig. [51 
In particular, our model recovers the dependence on the 



ratios {a, /3} and on the scale 9s and performs better than 
the other models studied in this paper. To simplifiy the fig- 
ures, we show in Fig. |S]only the second-best model, i.e. the 
"J2.nl" ansatz (l?7l) using the effective kernel F2 nl from 
iScoccimarro fc C^chmani (pOOl and the power from our 
model (as in the upper blue dashed line in Fig. Other 
models show similar behaviors to those found in Fig. [5] for 
single-scale moments (i.e., a lack of power on moderate and 
small angular scales). As for the second-order statistics, the 
underestimate of the weak-lensing signal by the simulations 
appears at a smaller angle 9s than in the single-scale case 
shown in Fig. [5] because the factors a and /3 are larger than 
unity and increase the sensitivity to larger scales at fixed 
9s- 

9. Conclusion 

We have investigated the performance of current theoretical 
modeling of the 3D matter density distribution with respect 
to weak-lensing statistics, focusing on configuration-space 
statistics, specifically the convergence two-point and three- 
point correlation functions and the second- and third-order 
moments of the smoothed convergence and of the aperture 
mass. As in paper I, where we studied Fourier-space statis- 
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Fig. 7. Two-scale second-order moments {Ks{ds)Ks{a9s)) 
(upper panel) and (Map(0s)-^ap(a^s)} (lower panel) as a 
fmiction oi Og, at Zs — 1. We show the cases a = 2, 5, and 
10 from top to bottom. The symbols are the same as in 

Fig.m 



tics, we found that a model in troduced in previous works 
() Valageas fc Nishimichil [2011 al ibi), which combines the (re- 
summed) one-loop perturbation theory with a halo model, 
fares better than some other recipes based on fitting formu- 
lae to numerical simulations or more phenomenological ap- 
proaches. It yields a reasonable agreement with numerical 
simulations and provides a competitive approach, because 
it remains difficult and time-consuming to describe a range 
of scales that spans three orders of magnitude or more by 
ray-tracing simulations. 

One advantage of our approach compared with numer- 
ical simulations or fitting formulas is that it allows us to 
decompose the integrated weak-lensing signal over several 
contributions that are associated with specific properties of 
the underlying 3D density field. Thus, we can distinguish 
perturbative terms, which can be derived from perturbation 
theory, from nonperturbative terms that are associated for 
instance with one-halo contributions, which depend on the 
density profile and mass function of virialized halos. This 
is useful because i) these different terms suffer from differ- 
ent theoretical uncertainties and ii) it allows one to under- 
stand which aspects of the matter distribution are probed 
by weak-lensing statistics, while angular scales vary. 

Like for the Fourier-space statistics studied in paper 
I and t he 3D statistics studied in Valageas & NishimichU 
(|2011allbh . we found that including one-loop terms in the 



Fig. 8. Three-scale third-order moments at = 1 as a 
function of ^s- Upper panel: {K.s{9s)ns{oi0s)ns{P0s))- Lower 
panel: {M^p{9s)Mi,p{a9s)Ms.p{(i9s))- We show the cases 
{a, /?} = {2, 5}, {3, 9}, and {5, 10} from top to bottom. The 
blue dashed line is the "-F2.NT," ansatz (l?7l). using t he ef- 
fective kernel i^2,NL from lScoccimarro fc CouchmanI (|200lD 
and the power from our model (as in the upper blue dashed 
line in Fig. [2]), while the red solid line is our model. 



perturbative contribution brings a more significant im- 
provement compared with the lowest-order perturbation 
theory for three-point statistics than for two-point statis- 
tics. Then, while large scales are described by these pertur- 
bative contributions and small scales by one-halo contri- 
butions, the nonperturbative two-halo term that gives an 
additional contribution to three-point statistics is always 
subdominant. This is a nice property because this mixed 
term is more difficult to model and may be less accurate 
than other contributions (see also paper I). 

Consequently, our model provides reliable predictions 
for weak-lensing statistics, from small to large scales, and 
for a variety of cosmologies. It could still be improved 
in various manners. First, the accuracy of the perturba- 
tive contribution may be increased by including higher or- 
ders beyond one-loop or by using alternative resummation 
schemes. Second, the underlyin g halo mode l could be re- 
fined to include substructures (Shcth 2003; Giocoli et alj 
. 2010.). deviations f rom spherical profiles (Jing fc Suto 200 2| 
Smith et"aI]|2006D . or the effect of baryons (jGuillet et al.l 
2010D . Next, the model could be generalized to non- 
Gaussian initial conditi ons, which yield disti nctive signa- 
tures in the bispectrum (jSefusatti et al.|[20Tol) . 
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